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Abstract 
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1 Introduction 


The Kibble-Zurek mechanism eng characterizes the spontaneous formation of defects in the 
process of system equilibration from an initial strongly nonequilibrium symmetric state to 
an equilibrium state with broken symmetry. This is a very general mechanism illustrated by 
numerous examples ranging from various condensed-matter systems to cosmological objects 
ME]. It has been experimentally observed, in ion crystals |6l[7|, superconducting dims [8], 
superfluid ^He and superfluid ^He, modeling cosmological string formation MI], in trapped 
atoms exhibiting Bose-Einstein condensation |l2l|13], and in a homogeneous Bose gas |14j . 

The case of nonequilibrium Bose-Einstein condensation in dilute Bose gases has attracted 
special attention since the work by Levich and Yakhot [12] followed by many other publica¬ 
tions studying, both theoretically [T6l - [T8] and in computer simulations [T9H2T], the process 
of equilibration of an uncondensed system to its Bose-condensed state. The general un¬ 
derstanding of the process of this equilibration, passing through the states with topological 
defects, is as follows. 

One starts with a strongly nonequilibrium state possessing no Bose-Einstein condensate, 
which is termed wave turbulence or weak turbulence. In the process of equilibration, the 
Kibble-Zurek mechanism is exhibited, hrst, by the appearance of local superfluid regions of 
condensate, with the characteristic size of coherence length, being yet not correlated with 
each other. These regions, called by Kibble BE] protodomains or cells, in the case of Bose 
condensation, are the germs, or grains, or droplets of Bose-Einstein condensed atoms inside 
uncondensed surrounding. Such a granular state represents the stage of strong turbulence, 
which can be called grain turbulence in order to stress the physical origin of the objects form¬ 
ing the state. This is similar to naming weak turbulence as wave turbulence for emphasizing 
that it is formed by small-amplitude waves. 

In the next stage of equilibration, the number of grains, or cells, grows so that they 
start fusing with each other forming numerous vortices [I1I21II2]. When no anisotropy is 
imposed on the system, the created vortices form a random tangle corresponding to vortex 
turbulence [22] . 

In the process of the following equilibration, the number of vortices diminishes, and the 
system passes to a vortex state with several vortices, whose number continues decaying, so 
that hnally the system relaxes to its equilibrium Bose-condensed state without vortices. 

In the present paper, we demonstrate, both experimentally and by computer simulations, 
that it is feasible to realize the inverse Kibble-Zurek scenario, starting with an equilibrium 
Bose-Einstein condensate, exciting which, we go through the sequence of states with sponta¬ 
neously arising topological defects, as in the equilibration process, but in the inverse order. 
Perturbing the system, we subsequently transfer it into the states of weak nonequilibrium, 
with no vortices, then to a vortex state, with several vortices, then to the states of vor¬ 
tex turbulence, grain turbulence, and, hnally, to the state of wave turbulence, where the 
condensate is completely destroyed. 

This is demonstrated by considering a trapped Bose-Einstein condensate, perturbing 
which we realize the whole sequence of states whose order is opposite to that arising in 
the process of equilibration of an uncondensed system to its Bose condensed state. All 
intermediate states with topological defects are recovered, although in the opposite to the 
Kibble-Zurek mechanism order. 

Although we call the process of transferring a Bose gas from its equilibrium state to a 
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strongly nonequilibrium state inverse Kibble-Zurek scenario, this does not necessarily imply 
that all particular details of this mechanism have to be exactly the same, but in the opposite 
order. Moreover, there is an important difference. One often considers the direct Kibble- 
Zurek mechanism caused by quenching the temperature of an external bath. However, in our 
case, there is no external bath hxing or varying temperature. We realize the transfer from an 
equilibrium state to highly nonequilibrium states by strong time-dependent perturbations 
of a finite closed system that is not connected to any bath. So, our system is strongly 
nonequilibrium. This principally distinguishes it from the situation of a stationary system, 
where a kind of vortices can be thermally excited in the vicinity of the Bose-condensation 
temperature [23H26]. These thermally excited vortices in a stationary system cannot form 
turbulence, while in our case of a strongly nonequilibrium system, vortex turbulence does 
happen. 

Naming the process, we describe, as the inverse Kibble-Zurek scenario, we aim at em¬ 
phasizing that the created nonequilibrium states with topological defects occur in the order 
opposite to the direct Kibble-Zurek mechanism. However, this does not mean the exact cor¬ 
respondence to all particularities of this mechanism. More precisely, the process, we realize, 
is inverse with respect to the so-called nonequilibrium Bose-Einstein condensation, because 
of which it could be termed the inverse Bose-Einstein condensation. However, we wish to 
stress the general possibility of reverting, not only Bose condensation, by also other different 
nonequilibrium phase transitions, realizing inverse nonequilibrium transitions. 


2 Perturbation of Bose-Einstein condensate 


We employ the same experimental setup as in the experiments studying quantum turbulence 
of trapped ®^Rb atoms [271 - 12^ . The gas of ®^Rb atoms is cooled down to the state where 
almost all atoms, up to 70% are condensed, with the number of atoms in Bose-Einstein 
condensate being N ^ 2 x 10^. The harmonic trap has cylindrical shape, with characteristic 
frequencies a;_L = 27r x 210 Hz and ooz = 2 tt x 23 Hz. Then the aspect ratio 


a = — - 
“ u± 



( 1 ) 


where l± and Iz are the transverse and axial oscillator lengths, is a = 0.11. 

After forming equilibrium Bose-Einstein condensate, the atomic cloud is subject to an 
oscillatory trap modulation that shakes the condensate without imposing a rotation axis, but 
preserving the cloud isotropy [SO]- As in our previous papers [271 - I2U] . we use the perturbing 
potential 


HI 

V{r,t) = -Ql{t){x'-x',r + 

in which m is atomic mass, x', y', z', x'q, y^, Zq are the variables shifted and tilted with respect 
to the real spatial variables x, y, z, and fla{t) = AaU)a[l — cos{ut)] are alternating frequencies 
oscillating with an effective amplitude and frequency a; = 27r x 200 Hz. This perturb¬ 
ing potential distorts the atomic cloud with a combination of deformations, rotations, and 
displacements [271 - I29] . The modulation amplitude and time can be varied in a wide range. 


m 




( 2 ) 
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thus, varying the energy injected into the system. The injected energy per atom can be 
evaluated [211132] as 



expressed through the energy Hamiltonian 

^ = j ^ dr + ^^0 J lh(r,t)|^ dr . (4) 

Here r] = 'r]{r,t) is the condensate wave function, U = U{r,t) is the total external potential 

U{r,t) = U{r) + V{r,t) , (5) 

including a trap potential U{r) and the modulation potential V (r, t). The atomic interaction 
strength is 

$0 = dvr — , (6) 

m 

with Qg being scattering length and m, atomic mass. We use the system of units, where the 
Planck constant h is set to one. 

The atomic cloud at initial time is assumed to be almost completely condensed, so that 
the condensate function is normalized to the total number of atoms 

= y \ri{r,0)\‘^ dr . (7) 

The equation of motion for the condensate wave function is 

' t ^ ^ 

Increasing the amount of the injected energy Einj disturbs the ground-state condensate 
and generates a sequence of nonequilibrium states. Although the created nonequilibrium 
states can be compared to some equilibrium states [33|, strictly speaking, they are quite dif¬ 
ferent. As is emphasized in Introduction, in equilibrium systems, there can be no turbulence, 
while in our hnite closed system turbulent motion does appear. 

Observations were done in time-of-flight experiments, after switching off the trapping 
potential. In the typical time-of-flight procedure, the chaotization [33] and nonadiabatic [33] 
effects do not occur and the expanding atomic cloud preserves the features of the state inside 
the trap [27H29]. 

We also accomplished numerical simulations for the whole process of producing strongly 
nonequilibrium states, with the setup and all parameters being exactly the same as in the 
experiment with ®^Rb. This was done by solving the nonlinear Schrodinger equation (8) by 
the methods of static [33] and dynamic [33] simulations. The dissipation is taken into account 
by introducing the phenomenological factor 7 [381139]. This dissipation slightly diminishes 
the number of atom in time, although not more than 10%. Increasing the amount of the 
injected energy, we have reproduced the whole sequence of nonequilibrium states: weak 
nonequilibrium, vortex state, vortex turbulence, grain turbulence, and wave turbulence. All 
numerical results are in good agreement with experimental data. 

By means of the alternating perturbation potential, modulating the trap, we inject into 
the trap the energy per atom Einj. By increasing the amount of the injected energy, we 
observe the following sequence of nonequilibrium states. 
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2.1 Weak nonequilibrium 

Under weak perturbations, only elementary collective excitations arise, with the cloud bend¬ 
ing and center-of-mass dipole oscillations, but without topological defects. Such elementary 
excitations are typical of collective excitations arising in weakly perturbed hnite quantum 
systems [31]. These excitations in Bose gases have been widely studied [3D|- So that we do 
not need to spend time for their discussion. 

In this regime, the average kinetic energy is yet much lower than the interaction energy. 
The effective interaction parameter is dehned as 

^ dvriV ^ . (9) 

This parameter is large because of the large number of atoms N, allowing for the use of the 
Thomas-Fermi approximation, where the effective atomic cloud radius and length are 

/15 \ 

R = i — ag j l± = 1.036 l±{agy^^ , aL = 2R = 2.072 l±{agY^^ . (10) 

We shall use these expressions in our following estimates. 

2.2 Vortex state 


The regime of weak perturbation lasts till the injected energy reaches the value 2Eyor 
required for generating a vortex-antivortex pair. Under the employed perturbation potential, 
imposing no rotation axis, vortices usually arise in pairs or larger groups, such that the total 
rotation moment be close to zero. 

The energy of a vortex, per atom, in the Thomas-Fermi approximation, m. reads as 

Eyor = 2mR^ (3.364 R^/^) , (11) 

where p is atomic density. Depending on the value of the density, the number under the 
logarithm can be slightly different, which, however, does not change much the vortex energy. 
For example, the average density can be dehned as N/ttR^L. But for a vortex along the axis 
of the trap, the central density is 

15iV ^ aN 

^ ^ ^ l\{agfR ■ 



When there are several vortices, the related densities p are different, because the density of 
the atomic cloud in a trap varies in space |40] - [42] . However, this variation does not essentially 
change the related vortex energy, since the density is inside the logarithm. Therefore, to 
estimate the vortex energy, we can accept density (12), which gives 


Eyor = 0.93 -—^ XniQAAag) . (13) 

[agyA 

In the present setup, a = 0.11 and g = 1.96 x 10^. The transverse trap energy E± = 
hjjj± = 0.869 X 10“^^ eV, which gives Eyor = 0.257 x 10“^^ eV. In units of the transverse trap 
energy, Ey^r = 0.296£'j_. Hence, 2Eyor = 0.591i?j_. 
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After the injected energy surpasses 2Eyor, the acting perturbation starts generating vor¬ 
tices. Since no single rotation axis is imposed onto the system, the creation of vortices and 
anti-vortices is equally likely. The typical situation is that vortices and anti-vortices appear 
together in groups [271j29] . The vortex-antivortex pairs do not form bound dipoles, being 
not correlated with each other, except having opposite vorticities. As we have checked by 
numerical simulations, practically all vortices have vorticity one or minus one, so that the 
total momentum is close to zero. The vorticity was calculated by the method described in 
Ref. [13] . 


2.3 Vortex turbulence 


Pumping into the trap more energy creates more and more vortices and anti-vortices. The 
used perturbing potential does not impose a rotation axis, because of which the increasing 
number of vortices does not lead to the creation of a vortex lattice, but results in the 
formation of a random vortex tangle representing, according to Feynman [22], quantum 
vortex turbulence [38] 13^. In this regime, the vorticity vectors are randomly oriented [44] . 
This regime develops, when the number of vortices V is such that the mean distance between 
their surfaces becomes twice the coherence length ^ ~ h/mc, where c ~ {h/m)^J4:'Kpas is 
sound velocity. This gives 

/ \ 2 


Nc 



(14) 


where r± is the atomic cloud transverse radius. For the treated case of ®^Rb atoms, ^ ~ 
2 X 10“® cm and r_L 4 x 10“^ cm. Then this estimate gives V ~ 25, which is in agreement 
with both numerical simulations and experiment. 

The number of vortices as a function of time, that is, as a function of the injected energy, 
under a fixed perturbation amplitude, is shown in Fig. 1, found in numerical simulations. As 
is seen, when the number of vortices reaches about 25 there occurs a change in the behavior 
of the function N{t), which is connected with the arising vortex turbulence. The following 
changes in the function N{t) will be discussed below. 

In the time-of-flight observation, the turbulent atomic cloud expands isotropically [27f[29] . 
This suggests that the observed random vortex turbulence is a kind of Vinen turbulence 
[38l[39|. Although a more precise classification of the turbulence type requires the study of 
spectra, which is yet in progress. 

The energy required for creating the random vortex turbulence can be estimated as 
Etur ~ NcEyor ~ 0.643 X 10“^^ eV, or in units of the trap energy, Etur = 7AE±. In the 
regime of vortex turbulence, the energy of a vortex is yet lower than the energy of vortex 
interactions. 


2.4 Grain turbulence 

When the number of vortices in the random tangle grows to the critical number N*, the 
vortex interactions become very important. Then vortices start strongly interacting and 
destroying each other. This happens when the interaction energy between two vortices, at 
a typical distance 6, reaches the energy of a vortex [381140] . This critical number of vortices 
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can be estimated as 

iv; = . (15) 

In the present setnp, the energy of a vortex becomes eqnal to the interaction energy of two 
vortices at the distance S ~ 0.71 x 10“^ cm, which yields N* ^ 30. This estimate is again 
in good agreement with nnmerical simnlations and experiment. 

The grains have the linear sizes between 10“^ cm and 5 x 10“® cm, which defines the 
typical size of the grain as 3 x 10“^ cm. This typical size is of order of the coherence length 
as it shonld be for a coherent droplet. Their density is mnch larger than that of their 
snrronnding. They look like the droplets of water flowing inside fog. 

The corresponding injected energy, needed for prodncing grain tnrbnlence, can be esti¬ 
mated as Efog = N*E^ori which eqnals Efog = 0.771 x 10“^^ eV, that is, Ej^g = S.S7E±. 

Fignre 2 presents the seqnence of experimentally realized states, observed in time-of-flight 
experiments with ®’^Rb: a state with several vortices, vortex tnrbnlent state, and the density 
confignration, typical of a granular state. Brighter colour corresponds to higher density. 

Numerical simulations confirm that the grains have the sizes close to the coherence length 
^ ~ 10“® cm and that the phase inside a grain is constant, hence the grains are really 
Bose-condensed coherent droplets. If the pumping by the alternating external potential is 
stopped, each grain survives during the lifetime 10“^ s, which is much longer than the local 
equilibration time of 10“^ s. That is, the fog made of these Bose-condensed grains can be 
treated as a metastable object, similarly to the fog made of water droplets [IS]. The density 
inside a grain is up to 100 times larger than that of their surrounding. It is worth stressing 
that the grains as well as their surrounding are formed by the same atoms. In that sense, 
the system should not be treated as composed of different species [16], but rather it could 
be interpreted as consisting of different phases [33] . 

During the regime of grain turbulence, the grains do not essentially change their sizes, 
increasing only very little. Since the transition between vortex turbulence and grain turbu¬ 
lence is a crossover, in the grain-turbulent state there can occur occasional vortices, although 
they are rare. 

2.5 Wave turbulence 

After the injected energy becomes comparable to the energy Ec = ksTc, where is the Bose- 
Einstein condensation temperature, all condensate is destroyed, and the system represents 
the state of wave turbulence, or weak turbulence. For ®^Rb atoms, Tc = 2.76 x 10“^ K. 
Hence, E^, = 0.238 x 10“^° eV, which gives Ec = 27AE±. 

Contrary to grains, the waves in the regime of wave turbulence are small-amplitude 
nonuniformities, being just three times denser than their surrounding. The typical size of a 
wave is of order 10“^ cm, and the phase inside it is random, which implies that there are no 
coherent regions that would be related to Bose-Einstein condensate. In the regime of wave 
turbulence, the Bose condensate is getting destroyed. 

Because of the high energy required for creating the wave turbulence, this regime has not 
yet been achieved in experiments, but has been realized in numerical modeling. The whole 
sequence of nonequilibrium states, found numerically, is presented in Fig. 3. 

Recall that the transitions between the dynamic regimes are crossovers, but not sharp 
phase transitions. But, as is seen from Fig. 1 for the number of vortices as a function of time. 
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the transitions from the vortex state to vortex tnrbulence and from vortex tnrbnlence to grain 
tnrbnlence are rather sharp. Contrary to this, the transformation from grain tnrbnlence to 
wave tnrbulence is rather gradual, the grains and waves coexist in a wide temporal range. 
The number of coherent grains diminishes with time, while the number of waves increases. 
The system coherence is being destroyed in a slow pace. First, coherence is destroyed at the 
surface of the atomic cloud. The snapshots of the phase are shown in Fig. 4. The crossover 
line, conditionally separating the regimes of grain turbulence and wave turbulence, can be 
located at the moment, when the system coherence, associated with grains, is destroyed in 
approximately half of the system, the incoherent region being associated with waves. 

During the process of pumping, the kinetic energy of the system, Ekin, increases, as 
compared to the energy of atomic interactions Eint- This is illustrated in Fig. 5. In the 
region of weak nonequilibrium, the interaction energy is much larger than kinetic energy. 
The vortex state, appearing at t^or = 18 ms and ending at ttur = 28 ms, is characterized by 
the ratio of kinetic to interaction energy between 1.8 and 4.2. Vortex turbulence, arising at 
ttur = 28 ms and lasting till tf^g = 37 ms, is described by the ratio E^in/E^nt between 4.2 
and 6.7. Under grain turbulence, existing in the interval from tjog = 37 ms to = 120 ms, 
the ratio Ekin/Etnt varies from 6.7 to 74. In the regime of wave turbulence, kinetic energy 
is much higher than interaction energy, so that Ekin/Eint is of order 50 to 100. 

3 Discussion 

When a system, being prepared in a strongly nonequilibrium symmetric state, equilibrates 
on a hnite time scale to a state with broken symmetry, there appear topological defects, 
such as vortices, protodomains, cells, or grains, strings, and monopoles. This is called the 
Kibble-Zurek mechanism that can be illustrated by various systems, from different kinds of 
condensed matter to cosmological objects [THSl ITTlITS] . In the particular case of nonequilib¬ 
rium Bose condensation, the equilibration process starts with the nonequilibrium state of 
wave turbulence and ends in an equilibrium Bose-condensed state [4^150] . 

We have shown that there exists the inverse Kibble-Zurek scenario, when we start with an 
equilibrium system with broken symmetry and, by imposing perturbations, transform it to a 
strongly nonequilibrium symmetric state through the sequence of states with spontaneously 
arising topological defects. This is demonstrated by considering a trapped Bose-Einstein 
condensate, perturbing which we realize a sequence of states whose order is opposite to that 
arising in the process of equilibration of an uncondensed system to its Bose condensed state. 
All intermediate states with topological defects are recovered, although in the opposite to the 
Kibble-Zurek mechanism order. Starting with an equilibrium Bose condensate, by increasing 
perturbation, we consecutively observe the creation of separate vortices, vortex turbulence, 
condensed cells, or grains, and wave turbulence with Bose condensate being destroyed. These 
are the same states that would appear in the opposite order, if an uncondensed nonequi- 
librium system would equilibrate to its Bose condensed equilibrium. We demonstrate the 
inverse Kibble-Zurek scenario both experimentally, by perturbing the Bose-Einstein conden¬ 
sate of trapped ®^Rb atoms, and also by accomplishing numerical simulations for the same 
setup as in the experiment, the experimental and numerical results being in good agreement 
with each other. 

The summary of the realized states, depending on the amount of the injected energy, is 


as follows: 


^inj 0 

0 < Einj < 2E vor 

^EuQP ^ Ej^Ylj ^ El^rj. 
Etur ^ Einj ^ ^fog 
^fog ^ ^inj ^ E(^ 
Einj E(^ 


(equilibrium state) 
(weak nonequilibrium) 
(vortex state) 

(vortex turbulence) 
(grain turbulence) 
(wave turbulence) . 


In terms of temporal scale, the observed nonequilibrium regimes occur in the intervals 


t = 0 

0 i tvOV 
ivor ^ tfnr 

ttur ^ t ^fog 
tfog ^ t <i tc 

t > tc 


(equilibrium state) 
(weak nonequilibrium) 
(vortex state) 

(vortex turbulence) 
(grain turbulence) 
(wave turbulence) , 


where Eor = 18 ms, ttur = 28 ms, tfog = 37 ms, and tc = 120 ms. 
Introducing the relative injected energy 






inj 


E„ 




normalized to the vortex energy, the sequence of the states in the inverse Kibble-Zurek 
scenario, we have realized, is characterized by the relations 


e = 0 
0 < e <2 
2<e <25 
25 < e < 30 
30 < e < 100 
e > 100 


(equilibrium state) 
(weak nonequilibrium) 
(vortex state) 

(vortex turbulence) 
(grain tnrbulence) 
(wave tnrbulence) . 


The injected energy (3), for the used oscillating perturbation potential is approximately 
proportional to time, after it is much longer than the period of the modulating held, t ^ 
27i/u, which is of order 10 ms. Thus, we have 


^fog 


E, 


tur 


^fog 


^tur 


1 , 


Ec tc 

- r\j - 

^fog tfog 


The different dynamic regimes are not separated by sharp bonndaries, bnt, strictly speak¬ 
ing, occnr continnously. That is, these dynamic transitions are crossovers, bnt not abrupt 
phase transitions, althongh they may happen rather noticeably. Therefore the neighboring 
states can occasionally display the features typical of their lower-energy member. For in¬ 
stance, in the regime of grain turbulence, there may happen occasional vortices. And in the 
regime of wave tnrbulence, waves coexist with grains. 

Thus we have demonstrated that the Kibble-Zurek mechanism is reversible, at least 
for the process of Bose-Einstein condensation versus Bose condensate destrnction. The 
demonstration of the inverse Kibble-Zurek scenario for trapped atoms opens a new direction 
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of studying the universality of the reversibility of this mechanism, investigating whether it 
is reversible for other physical systems of different nature. 
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Figure Captions 


Figure 1. Number of vortices as a function of the pumping time, under fixed modulation 
amplitude, found in numerical simulations. 


Figure 2. The sequence of nonequilibrium states observed in experiment: (a) vortex 
state, (b) vortex turbulence, (c) grain turbulence. The density snapshots, obtained in the 
time-of-flight setup, are presented. Darker spots correspond to lower density. 


Figure 3. The sequence of nonequilibrium states realized in numerical modeling: (a) 
vortex state, (b) vortex turbulence, (c) grain turbulence, (d) wave turbulence. The density 
cross-sections are demonstrated. Brighter colour corresponds to higher density. 


Figure 4. The sequence of snapshots of the system phase at consecutive moments of 
time, the phase being gradually destroyed under increasing injected energy. 


Figure 5. The ratios of interaction energy to kinetic energy of the system, Eint/Ekin 
(thick line), and of kinetic to interaction energy, E^in/Eint (thin line), as functions of the 
time of pumping. 


14 




Figure 2: The sequence of nonequilibrium states observed in experiment: (a) vortex state, 
(b) vortex turbulence, (c) grain turbulence. The density snapshots, obtained in the time-of- 
flight setup, are presented. Darker spots correspond to lower density. 
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Figure 4: The sequence of snapshots of the system phase at consecutive moments of time, 
the phase being gradually destroyed under increasing injected energy. 


Figure 3: The sequence of nonequilibrium states realized in numerical modeling: (a) vortex 
state, (b) vortex turbulence, (c) grain turbulence, (d) wave turbulence. The density cross- 
sections are demonstrated. Brighter colour corresponds to higher density. 
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Figure 5: The ratios of interaction energy to kinetic energy of the system, Eint/Ekin (thick 
line), and of kinetic to interaction energy, Ekin/Eint (thin line), as functions of the time of 
pumping. 
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